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I. INTRODUCTION 



Symmetry reduced models of General Relativity have been extensively studied and ap- 
plied in the framework of cosmology In cosmological models the infinite degrees of freedom 
of the full theory are reduced to a finite number describing the large scale structure of the 
universe. Current observations seem to be best described by (spatially flat) homogeneous 
and isotropic Friedman-Robertson- Walker space-times. For that reason, the assumption of 
homogeneity and isotropy is by now an accepted basic ingredient in the so-called standard 
model of cosmology. 

Under quite general conditions, general relativity predicts the development of space-time 
singularities. One prime example of such singular behavior is the presence of the initial Big- 
Bang singularity in cosmological models. Energy densities as well as space time curvature 
diverge approaching the initial Big-Bang singularity, and any description based on classical 
space time notions becomes inapplicable. General Relativity is no longer valid and has to 
be replaced by a quantum theory of gravity. 

Although no complete quantum theory of gravity is yet available, loop quantum gravity 
appears as an promising candidate [l|, Q, 0, Qj- Loop quantum gravity is a background 
independent non perturbative approach to the quantization of general relativity. The main 
prediction of the theory is the discreteness of geometry at the fundamental level: operators 
corresponding to geometric quantities such as the area and volume have discrete spectra 



HQ. 



Geometry is quantized and at small scales the smooth notion of space and time 
simply can no longer be applied. One of the main clear-cut applications of the theory, where 
the discovery of discreteness plays an important role, is the description of the fundamental 
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131 ] - These fundamental 



degrees of freedom of black hole (isolated) horizons [ 
excitations give rise to the correct value of black hole entropy that is expected from semi- 
classical considerations. 

The application of loop quantum gravity to cosmological models is known as loop quan- 
tum cosmology (LQC) 14] (See also |l^ for an alternative approach). The simplification 
arising from the symmetry reduction makes it possible to address physical questions that 
still remain open in the full theory. For instance LQC provides a novel paradigm for the 
understanding of the fate of the Big-Bang singularity in quantum gravity. The important 
physical question of whether classical singularities of general relativity are not present in 
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the quantum theory is answered in the affirmative in LQC 16]. The basic mechanisms lead- 
ing to this striking result can be traced back to the fundamental discreteness of geometric 
operators which drastically changes the quantum evolution equation in the deep Planckian 
regime. In addition, consistency conditions,known as the dynamical initial condition, seem 
to severely restrict the freedom in the choice of 'initial' values for the evolution and it is 
argued that these conditions may single out a unique wave function of the universe [1 71]. It 
has also been shown that loop quantum cosmology leads to modifications of the classical 
Friedman equation. Recent work has shown that these effects can provide the proper initial 
conditions for chaotic inflation and may lead to measurable effects on the CMB [18]. 

In spite of these successes, a complete characterization of the physical Hilbert space of 
loop quantum cosmology remains an open issue 1 . The definition of the physical inner product 
is the missing piece of the puzzle in LQC. Without the physical inner product, there is no 
notion of a probability measure which is an essential ingredient for physical prediction in 
the quantum theory. Any interpretation of the wave function as a probability amplitude is 
futile without the probability measure provided by the physical inner product. 

In addition the physical inner product provides a means for eliminating pathological 
solutions to the constraints. For instance it can be shown that the space of solutions of 
LQC is generally non separable, i.e., an uncountable infinity of solutions to the quantum 
constraint equation exist. Such a huge physical Hilbert space is certainly not expected for 
a system with finite number of degrees of freedom. Some of these solutions exhibit peculiar 
properties such as wild oscillatory or diverging behavior in the regime where one would 
expect agreement with standard Wheeler-DeWitt quantum cosmology. In the absence of 
the physical inner product, so far spurious solutions are eliminated using the dynamical 
initial condition coupled with heuristic arguments based on semi-classicality requirements. 
Although these requirements are physically motivated, unphysical solutions can be easily 
identified if a notion of physical inner product is provided. More precisely, on the basis 
of the notion of physical probability physical states are defined by equivalence classes of 
solutions up to zero norm states. It is hoped that the ill-behaved solutions mentioned 
above will be factored out by this means once a suitable notion of physical inner product is 

1 This remains an open issue in the full theory of LQG. The master constraint approach |l9j and the spin 
foam formulation [20j are current proposals to solve this problem. 
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provided. 

There exists a simple method known as group averaging to provide a definition of the 
physical inner product in generally covariant systems when the constraints generate a group 



action 21]. In this paper we review and use the technique called refined algebraic quanti- 
zation 22j as a recipe to construct the physical Hilbert space. The physical inner product 
is constructed in a simple manner if one can define a unitary operator generating 'transla- 
tions' along the gauge orbits generated by the quantum constraints in a suitable kinematical 
Hilbert space. This technique cannot be applied to the standard formulations of loop quan- 
tum cosmology because the quantum Hamiltonian constraint is not a self adjoint operator 
(in the kinematical Hilbert space) and thus it does not lead to any well defined unitary 
'evolution' when exponentiated. 

In this paper we present a model of loop quantum gravity that arises from a symmetry 
reduction of the self dual Plebanski action. In this formulation the symmetry reduction leads 
to a very simple Hamiltonian constraint that can be quantized explicitly in the framework of 
loop quantum cosmology. The model is defined in the Riemannian sector which constitutes 
a limitation when considering physical predictions. However, we will present a method for 
transforming to the Lorentzian sector and then argue that the main results of the model 
are manifested in the Lorentzian sector. The model shares the same kinematical Hilbert 
space with standard loop quantum cosmology, but differs in that the resulting Hamiltonian 
constraint is self-adjoint and thus the group averaging technique can be used to construct the 
physical inner product. When matter is incorporated in the form of a cosmological constant, 
we can solve the model exactly both for the physical wave functions and the physical inner 
product. The physical inner product reduces the infinite set of solutions of the quantum 
Hamiltonian constraint equation to one physically normalizable solution. This constitutes 
an interesting example where the very large (non separable) kinematical Hilbert space is 
reduced to a (trivially) separable physical Hilbert space. 

We end with a discussion of the relevance of the model to standard LQC We show that 
despite the simpler constraint the model shares the same qualitative features of standard 
LQC with one important difference. We will then discuss the implications that can be 
derived about LQC in general. 
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II. CLASSICAL SYMMETRY REDUCTION 



Our starting point is the classical self dual Riemannian gravitational action which, as 
noted by Plebanski[23], can be written as a constrained SU(2) BF theory 

1 
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AS, AS, (1) 



where £ = £* dx^dx u Ti is an SU(2) Lie algebra valued two form, A = A^dx^Ti is an SU(2) 
Lie algebra valued connection, r, are the generators of the Lie algebra of SU(2), and A is 
a cosmological constant. In this paper we set the Planck length l p = ^/8nGh equal to one 
and restrict the model to flat k=0 cosmologies. The tensor ipij is symmetric and acts as a 
Lagrange multiplier enforcing the constraint S* A S J — |<5 lJ E fc A E fc = 0. Once this constraint 
is solved, the action becomes equivalent to the self-dual action of general relativity. 

In order to reduce the action ((TJ to spatial homogeneity and isotropy, we write the action 
explicitly in terms of coordinates separating space and time as 

S GR [A,E,B] = Jdt J dh A\E* + A\D a E« + e ahc B\Fl c 



'0. 



'j 



-BlE aj + -B j n E ai - - 
2 a 2 a 3 



+ A5j£f (2) 



where we have introduced the definitions E ai = 2e abc T, l bc , B\ = 2£Q a , and e abc = e 0abc . 



Spatially homogeneous and isotropic connections can be characterized as| 

K = AkWa (3) 

where the parameter A encodes the gauge invariant part of the connection, A} is an SO (3) 
matrix, and uo 1 = u^dx a (7 = 1,2,3) are left invariant one forms with respect to the transla- 
tional symmetry associated with homogeneity. Since we are restricting ourselves to spatially 
flat cosmologies, the translation symmetry group is isomorphic to K 3 . Due to isotropy the 
time component of the connection A l Q vanishes identically 2 . The matrix A} is pure gauge 



2 According to the general analysis of symmetric connections for arbitrary symmetry groups |2'1| the time 
component of the connection is a connection on a symmetry reduced principal fiber bundle. This fiber 
bundle has as its gauge group the centralizer Z\ := Zg(X(F)) where G is the gauge group of the theory, 
F is the isotropy group, and A is a homomorphism from F to G. For our model the gauge group as well 
as the isotropy group is SU(2). The homomorphism A is thus the identity map and the centralizer Z\ 
only contains the identity. The time component of the connection thus vanishes. 
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and satisfies 



A)A< = 5] 

e ijk A\A 3 jA k K = e IJK . (4) 

Using left invariant vector fields ujj = ujjd a dual to u 1 we can write the momentum canoni- 
cally conjugate to the connection as 

E? = yfqoEAjujj (5) 

where the density weight of E? is provided by the left invariant metric q 0ab = LO a u^8ij. Once 
again, the gauge invariant part of the canonical momentum is given by a single parameter 
E. Similarly B l a is given by 

B\ = BA\u{. (6) 

Written in terms of the reduced variables the homogeneous and isotropic gravitational 
action becomes 

S GR [A, E, B] = Jdt Jd 3 x 3^/q^EA + 6y^BA 2 + 3^/q^ABE 

= V Jdt 3EA + 6BA 2 + 3 ABE (7) 



where Vq = jd?x^/%. The constraint term -0, 



B\E a i + B{E ai - \5 lj B k a El 



in the full action 



(J2J vanishes identically for the isotropic model which can be checked by a direct calculation. 
Technically Vq diverges, since the left invariant metric is constant and the integral is over a 
non compact manifold. To overcome this, we restrict the integral to a finite cell with volume 
Vq and absorb this factor into the variables as E — > V 2 ^ 3 E } A — > Vq^ 3 A, and B — > V^ 3 B, 
whence the action becomes 

S GR (A,E,B) = Jdt3EA + B(6A 2 + 3AE). (8) 

It is clear from the action (JHJ) that E is canonically conjugate to A with Poisson bracket 
{A, E} = |, and B acts as a Lagrange multiplier enforcing the Hamiltonian constraint 
H = which here for the gravitational part is just H = 6A 2 + 3AE. Counting degrees of 
freedom we have two dynamical variables A and E with one first class constraint H thus 
giving zero degrees of freedom as expected for isotropic general relativity with no matter. 
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We now show that the action (JHJ) is equivalent to the standard isotropic action written 
in terms of a single scale factor a and lapse N with metric ds 2 = N 2 dt 2 + a 2 dx 2 . In the full 
BF theory, the densitized metric is given by the Urbantke formula [3| 

y/99r = le ijk e^% a E^Et. (9) 
Reducing this to isotropy gives 

2B 2 

9oo = 



E 

9aa = y (10) 

We thus have the relationships between the two sets of variables 

E = 2 a 2 

B = Na. (11) 

In the full theory, the connection is a sum of the spin connection plus the extrinsic curvature. 
Since we are only considering spatially flat models the spin connection vanishes and the 
connection is equal to the extrinsic curvature which implies for the isotropic model 

A= — . (12) 
2N v ; 

With the conventions chosen the action (EI) becomes 



S, 



GR 



[dt — 


2 •■ 1 .9 

a a H — aa 


1 N 


2 



+ 6ANa 3 (13) 



which coincides with the standard isotropic action of general relativity up to a total deriva- 
tive. 



III. QUANTUM THEORY 

We wish to quantize the symmetry reduced theory in the same manner as loop quantum 
cosmology by using techniques from the full theory adapted to the symmetry. Thus, the 
Hamiltonian constraint is represented as an operator on a kinematical Hilbert space using 
holonomies and fluxes as the basic variables. Physical wave functions are those annihilated 
by the constraint operator. In addition refined algebraic quantization provides a means for 
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constructing the physical Hilbert space. In this section we start with a formal introduction 
to refined algebraic quantization (or the group averaging technique) used to calculate the 
physical inner product. Next we carry out the program for the cosmo logical model considered 
here. 

A. Refined Algebraic Quantization 

The goal of refined algebraic quantization is to find a method for describing the physical 
Hilbert space, the space of wave functions annihilated by the constraint and normalizable 
with a physical inner product. Refined algebraic quantization provides a means for deter- 
mining both. We start with the formal definitions and notation. The discussion will not 
serve as a general review of the method, but will instead be tailored to the single constraint 
system at hand. 

The idea of refined algebraic quantization is to build the physical Hilbert space TL p h ys 
starting from a kinematical Hilbert space TLkin acquired by quantizing the theory first ignor- 
ing the constraint. The constraint is represented as a self-adjoint operator on this kinemat- 
ical Hilbert space. As in standard Dirac quantization, the physical states are annihilated 
by the constraint operator. If the eigenstates of the constraint operator are not normaliz- 
able then the physical states lie outside of the kinematical Hilbert space. In this scenario 
refined algebraic quantization provides a recipe for constructing the physical inner product 
and characterizing the physical Hilbert space. 

For the case where the eigenstates of the constraint operator are not normalizable the 
physical inner product is constructed as follows. A dense subspace $ C Tikin of the kine- 
matical Hilbert space is chosen and solutions to the constraint equation then lie in the 
topological dual $* of $. We denote the action of on <3> using Dirac notation; namely, 
given an element G and an element \<j>) G $ the action is denoted (ip | 0). For the 
construction of the physical inner product, an anti-linear map P : $ — > $* is required such 
that given a state \ip ) £ (-P(V'o)l G is a solution to the constraint equation in the sense 
that (P(ipo) \H\4>o) = for any |</> ) G <E>. Thus P maps onto the kernel of the constraint 
operator which consists of elements of $*. Technically P is not a projector since P 2 is ill 
defined. The spaces involved satisfy the triple relation $ C Tikin C $*. 

With the map P we can now define the physical inner product and characterize the 
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physical Hilbert space. The physical inner product between two states {P((f>o) | , (P(ipo) I £ 
is given by 

(P(0 o )|P(^o)W:= Wo)]0o)- (14) 

where the l.h.s is just notation, and due to the anti-linearity of P the order has been reversed 
on the r.h.s.. A unique physical state is labeled by an equivalence class of states in $ C H-kin- 
Two states |-0 O ) and \ip' Q ) in <£> are equivalent if 

(Vol = (Vol + (^o| (15) 

for some \x ) G $ satisfying 

(P(x ) | x ) = 0, (16) 

ie. when (x | has zero physical norm. Under these conditions the physical inner product is 
independent of the state |Vo) £ $ used to represent the physical state in (P(Vo)| G The 
elements of the physical Hilbert space can therefore be labeled by an equivalence class of 
states in $ as defined above. 

So far, to construct the physical Hilbert space we have demanded the existence of the 
map P without providing a method of calculating it. We now describe the group averaging 
technique which provides such a method. The map is given by 



(P(0o)| = JdT{Ue~ lHT . (17) 

Clearly P maps onto the kernel of the Hamiltonian constraint, that is 

(P(Vo)|^|0o) = O (18) 

for any state \<f> ) G $. This can be shown by inserting a complete set of eigenstates of the 
constraint operator H\E n ) = E n \E n ) to get 

(P(Vo) \H\(f) ) = JdT (Vo | e~ ltlT H | O > 

= J2[dT(^ \ E n ) (E n \ e~ iitT H\<f>o) 

n J 

= $^(Vo| E n ) j dT e~ lE - T E n (E n | O > 

= ^(Vo I E n ) E n 5{E n ) (E n | O ) = (19) 
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where in the case where the eigenstates are not normalizable the sum over n is replaced by 
an integral. 

If the eigenstates are normalizable then the construction of the physical Hilbert space is 
straightforward. The map P now is a bona fide operator on the kinematical Hilbert space 
and behaves as a projector satisfying P 2 = P projecting onto the kernel of the constraint. 
The physical inner product becomes 

(P[0o]|P[V>o]W = (A\P\<f>o) (20) 
= (V>o I P 2 I 0o> (21) 

thus the physical inner product is equivalent to the kinematical one in the restriction of 
Tikin to the kernel of H . The definition of the equivalence relation remains the same with a 
simplification arising since P 2 = P. 



B. Kinematical Hilbert Space of Loop Quantum Cosmology 

We now turn to the cosmological model at hand. The first step required in refined 
algebraic quantization is to build the kinematical Hilbert space by ignoring any constraints. 
The kinematical Hilbert space for isotropic loop quantum cosmology has been constructed 
in [2^. We review the results here. 

Using ideas from the full theory, configuration variables are constructed from holonomies. 
For homogeneity and isotropy the holonomy algebra consists of the set of almost periodic 
functions; namely, those that can be written 

f(A) = h e 1 ^ (22) 
3 

where the sum contains a finite number of terms (j — 1, 2, ■ • • N and N < 00), fj £ C, and 
Uj £ R. The momentum variables consist of fluxes of the triad operator Ef on a 2-surface 
which after symmetry reduction are trivially proportional to the parameter E. The kinemat- 
ical Hilbert space is defined by representing this algebra using the Bohr compactification of 
the real line RboHt- RBohr is a compact Abelian group and is equipped with a Haar measure 
denoted by dfi with elements being labeled by A £ R. The Haar measure can be explicitly 
written as 

dfi= lim — / dA. (23) 
10 



The kinematical Hilbert space consists of functions which are square integrable with the Haar 
measure dfi which corresponds to the completion of the set of almost periodic functions in 
the norm of this measure. In the notation of section Ull Al $ is the space of almost periodic 
functions. With this measure an orthonormal basis consists of states \v) given by 

{A\v) = e iAv/2 (24) 

where the parameter v spans the entire real line. The basis states are normalizable in 
contrast to a standard quantum mechanical representation and satisfy the orthonormality 
condition 

(1/ | v) = <W (25) 

which can be shown using the measure ([23)1 . The basis states \v) are eigenstates of the triad 
operator E = —i\-§x'- 

E\v) = V - \v). (26) 

Geometrical operators can be built from the triad operator. In particular the volume oper- 
ator is given by V = \E\ 3 ^ 2 thus \u) represents a quantum state with definite volume given 
by (^r-) 3 / 2 . This provides the physical interpretation of the label v. Finally we can write the 
decomposition of the identity for both connection and triad representation as 



1 = J dfi \A) (A\ 

oo 

1 = £ (27) 

v=— oo 

where the sum over v is a continuous sum over all values of v on the real line. It is important 
to note that the kinematical Hilbert space here is non-separable being spanned by a non- 
countable set of basis states \v). 



C. Quantization of the Hamiltonian Constraint 

The homogeneous and isotropic model consists of the Hamiltonian constraint H c i ass = 
6A 2 + 3AE + -mH matter 3 . In representing the constraint as a self-adjoint operator on the 



3 The factor of 1/ \[E is added such that H ma tter agrees with the standard form of the matter Hamiltonian 
used in other formulations. Note that the constraint used here does not have units of energy hence the 
need for the extra factor. 
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kinematical Hilbert space no operator corresponding to the connection A exists; instead con- 
nection operators need to be represented in terms of the basic variables which are holonomies 
along edges. In particular the curvature term in the full action (J2J , which contributes the 
A 2 term after symmetry reduction, can be represented with holonomies around closed loops. 
Because of isotropy we can consider holonomies around squares on the manifold with edge 

1/3 

lengths chosen to be v V for some positive parameter is . The holonomy along an edge 
generated by the left invariant vector field uj is given by 

hj = exp(z/ AA 7 ) = cos(u A/2) + 2A/ sin(z/ A/2). (28) 

The 6BA 2 term in the action is then represented as 

6BA 2 -> -\ e w tr [Bkihjh K {hj)-\hK)~ x ] = ^ sm 2 (v A/2) cos 2 (z/ A/2). (29) 

In the limit where the closed loops are shrunk to a point by taking u to zero, the classical 
expression 6BA 2 is recovered. In the quantum theory since there is no operator correspond- 
ing to the connection A we cannot remove v§ and thus the parameter remains as a quantum 
ambiguity. In this paper we will not fix it's value, but it has been argued that it's value can 
be determined from the full theory of loop quantum gravity to be equal to v^3/4 j^J. 
The gravitational part of the Hamiltonian constraint thus is given by 

24 * 
H GR = — sin 2 (z/ A/2) cos 2 {u A/2) + 3AE . (30) 

To determine it's action on the basis states \v) we first consider the action of sin(z/ A/2) 
and cos(z/ ^4/2). Since we have (A | v) — exp(«74i//2) we have 



sin(z/ i/2) \v) = ^7 ( \v+J4>) - \v-Vq) 
cos(z/ i/2) \u) = -( \v+v ) + W-v ) ) (31) 
and thus the action of Hqr is 

H G r \v) = W+^o) - 2\u) + \u-4u ) ) + A^\u) . (32) 

The action of the curvature part of the constraint operator is thus seen to be a discrete 
approximation to an operator corresponding to the second derivative with respect to v and 
in the limit of small Uo approaches what would be the standard quantum mechanical operator 
6A 2 = -24|i 

12 



D. Path Integral Representation of P 



Now we explicitly compute the group averaging formula (JT7j) and give a definition of 
the map P. The matrix elements of the map P(v" , v') = iy" \P\v') — J dT (i>" \ e~ lHT \ v') 
can be written in terms of a path integral amplitude. The integrand looks like an ordinary 
quantum mechanical propagator and thus can discretized using the standard derivation. 
Decompositions of the identity ()27j) are inserted at n discrete time slices = ke with time 
step e = to get the propagator 

{u"\e- ikT \u') = limnjx; f^k e ^lL\MA k+1 -A k )/2-eH(A k ,u k ) 

fe=l l^fe^O 

= "J VAVEe iS " (33) 

with V\ = u' and u n+ x = u" and S being the discretized classical action S = Ylk=i ^(^fc+i - 
Ak)/2 — eH{Ak,Vk). Notice that the functional integration over the triad E is discretized 
as continuous sums over it's eigenvalues at the time slices Uk/6 and that H(Ak,iSk) = 
^§ sin 2 (z/ yW2) cos 2 (z/ A fc /2) + Au k /2 H — -^=H matter . We note that this path integral am- 
plitude is equivalent to the path integral amplitude arrived at after gauge fixing due to 



reparameterization invariance of the classical action 4 . 

Putting these ideas together, the matrix elements of the map P are calculated through 
the path integral as 

P(v", v 1 ) = [dT limn | E fak e iS \ . (34) 

fc=l (.v k ,vo ) 

We then use this to map any state l^o) = ^v^OvW) £ $ onto a solution of the constraint 
equation as in 



(PM\ = dT(ij \e 



-iHT 



Halliwell [2^ has considered the path integral for actions that are reparameterization invariant such as 
the action presented here. Starting with the path integral P = J VAVEVB exp[i f., dt3EA + BH] 
the one form component B is gauge fixed to be constant in time and after including ghost terms to 
make the path integral independent of the gauge choice the path integral becomes P — J dB(t" — 
t') JVAVEexp[i J*" dPZEA+BH] which with the redefinition T = B{t" -t') and t = B(t-t') takes on the 
form equivalent to the group averaging one P = J dT J VA VE exp[i J^dt 3EA+H] = J dT (v" \ e' iAT \ v') 
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u",u' 

= i^P(v",v'){v'\ ■ (35) 

u",u' 

The physical inner product between two physical states (P(0 O )| and (P(ip )\ is given by 

= m> i v") v i p i v) w 1 O ) 



P(z/ , Z/ ) 0o^' • 



(36) 



With the equivalence relations of equations (|T3jl and (fTB^) . we have a concrete recipe for 
characterizing the physical Hilbert space of the model. 

IV. EXAMPLE - COSMOLOGICAL CONSTANT 
A. Physical Hilbert Space 

As an example we consider the simplest non-trivial model with a positive cosmological 
constant and no other forms of matter. In this section we will show that this model can be 
solved exactly in the connection representation and we will explicitly construct the physical 
Hilbert space. In the following section we transform to the triad representation in order to 
discuss the relationship to loop quantum cosmology. 

We have already derived the form of the constraint which is given by 

24 



H GR = sinVo^/2) cosVo-4/2) + 3A£ 



24 



2 sin 2 (z/ A/2) cos 2 (u A/2) - iA 



_d_ 

dA 



(37) 



Because of the simplicity of the constraint operator we can solve exactly for the eigenstates 
which are given by solutions \E n ) satisfying H\E n ) = E n \E n ). This differential equation can 
be solved to give 



E n (A) = exp 



i j I E u 



cxp 



3 sin(2t/ -4) 



2kvl 



(38) 



where the eigenvalues E n are valued on the real line. The eigenstates are pure phase, but 
they are normalizable in the kinematical inner product. This is despite the fact that any 
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E n e K lies in the spectrum of ff. Since the eigenstates are normalizable we expect P to be 
a bona fide projection operator. The kinematical norm is 




and thus all of the eigenstates have kinematical norm equal to one. Following the general 
discussion of refined algebraic quantization, the physical Hilbert space consists of the space 
of zero 'energy' eigenstates and the physical inner product is equal to the kinematical one 
restricted to this space. We see that this space is one dimensional spanned by the solution 
\E Q ) given by 



We can also write out the matrix elements of the map P(A", A') = (A" \P\A'). Since the 
map projects onto the zero 'energy' eigenstate and satisfies P 2 = P the matrix elements 
must be of the form 



We thus see that despite the fact that the kinematical Hilbert space is very large, the 
physical Hilbert space is one- dimensional which is what we expect from a theory with zero 
degrees of freedom. In the next subsection we will see that in the triad representation there 
exist a non-countable set of solutions to the constraint equation — this is in fact generally 
the case in standard loop quantum cosmology — but we will show that the extra solutions 
all have zero physical norm and thus the physical Hilbert space is one-dimensional. 

B. Triad Representation 

In standard loop quantum cosmology a simple formula for the Hamiltonian constraint 
in the connection representation does not exist hence the constraint and solutions are more 
readily given in the triad representation. The space of solutions in standard loop quantum 
cosmology is non-separable whereas the physical Hilbert space of the model presented here 
is separable consisting of a single solution. Our aim in this section is to reconstruct the 
physical Hilbert space in the triad representation to facilitate a comparison with standard 
loop quantum cosmology. In particular the map P is calculated explicitly using the path 
integral and picks out a single state among all the solutions to the constraint equation. 




(40) 



P(A",A') = E (A")E (A') 



(41) 



15 



In the triad representation the action of the Hamiltonian constraint operator was given 
in equation (J32*|) 

H\u) = ( k+4z/ ) -2\u) + \u-4vq) ) + -v \v) . (42) 

We decompose a state lying in the dual of the kinematical Hilbert space as 

M = X>W- (43) 
The constraint equation W = leads to a difference equation for the coefficients ip u 

1pu+4u ~ + 1pu-4u = ^Az/q V 1p u . (44) 

It is immediately seen that the difference equation splits into an infinite number of isolated 
sectors with values of v separated by an integer times 4i/o. We label each sector by the 
parameter 5 G [0,4z/ ). Thus the sector labeled by 5 corresponds to values of v equal to 
. . . 5, 4:^0 + 5, 8^0 + 6 .... The difference equation for each sector is then a second order 
equation which can be solved exactly and the general solutions are given by 

3 \ ^ 5 , „ I 3 



- C x J u/{ 4u )+3/(2Aui) [2X^3 J +C 2 Y u/(4u )+3/(2Au3) [^X^ ) ^ 

where J and Y are the Bessel functions of the first and second kind respectively and C( 2 are 
constants. There are thus an infinite number of solutions determined by the infinite number 
of parameters Cf 2 for 5 G [0,4z/o). 

The matrix elements of the map P can be solved exactly for this model using the path 
integral (J3*lj) . The discretized action is 

T/e 

S = ^^(A k+1 -A k ) + e[6sm 2 {v A k )/it + Av k /2] 

k=0 

= i A k+t~ A k + eA] - 6esm 2 {u A k )/u 2 (46) 

fe=0 

thus the sums over v k can be performed since they are all of the form 
exp [iuk (Ak+i — Ak + eA) /2] which give delta functions. The result is 

„ .fr, , f°° dA k 5(A k+1 -A k + eA)e^" A -/ 2 e-^ A ^ 2 } 

(v"\e- lHT \v') = \imT] I J -°° V + 1 L (47) 

6-0ll\ e i£LU6"m>oA)A, 2 J 
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Let us then integrate over the Ak using the delta functions except for Aq and A n . This 
leaves us with one remaining delta function 

_ i& J^odAodA* 5(A n -A + TA)e-" A ^e-'^ 

exp 

The path integral is then given by 



T7> 



■M (rp _ sin(7V A) cos(2i/ A'-TVpA) 
o \ ^oA 



— oo 
oo 



3i 

7T? 



„, sin(7V A) cos^qA'-IVoA) 



P(z/, */) = / dA cM n e iu " An/2 e- iu ' Ao/2 dT 5(A"-A' + TA) 

) / 3isin(2yA") 3iA" \ ( 3isin(2v A') ZiA 1 

dA dA n e iu " An/2 e - iu ' Ao/2 e \ 2 "o A ^ ) e \ 2v Z A ^ 

— oo 

/oo 
dA Q dA n e iu " A " /2 e - iu ' Ao/2 E Q (A") E^A 7 ) (49) 
-oo 

which is just the Fourier transform of the map from the previous section. P thus maps onto 
the single normalizable zero 'energy' eigenstate \E ). 

We now formulate the physical Hilbert space by using the map (|49j) . P maps onto a single 
state in the triad representation which is given by the Fourier transform of the zero 'energy' 
eigenstate |-Eb)- Denoting \E ) = Yl v Pv\ v ) a calculation gives for the Fourier transform 

p = J J *V(M))+3/(2A^) (aJJjr) 4^ + 2X^ e Z ^ 



3 



i z. 



4u 

Remarkably, the path integral solution is non zero in precisely one sector 5 C corresponding to 
values of v where the Bessel function order ^ + is an integer. The value of the special 
sector is thus given by 5 C = 4z/ [1 — mod(^^, 1)] where mod(a, b) is the remainder of a/b. 
That the Fourier transform picks out certain modes can be seen from the fact that the state 
Eq(A) consists of a phase factor multiplied by a periodic function. The matrix elements of 
P in the triad representation are thus given by P vv < — (v \ P \ v') — p u jh7. Using this form of 
the map, any kinematical state will get projected onto the zero 'energy' eigenstate given in 
equation ([50)1 . The uncountable extra solutions (|4"5j) to the constraint equation in the triad 
representation are to be mod out of the physical Hilbert space. In this way the projector 
has reduced what was initially a non-separable space of solutions to the constraint equation 
to a one dimensional Hilbert space. This is because the generic solutions of the constraint 
correspond to zero physical norm states and therefore are equivalent to the zero state in the 
physical Hilbert space. 
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The previous statement requires a more precise explanation as the generic solutions (|45|) 
are not normalizable in terms of the measure (|23|). These solutions have zero physical norm 
in the following sense: The generic solutions (|45|) can be thought of as elements of $* defined 
in Section All Al Given a generic solution ip s G $* (as in (J45|) ) with values only in the sector 
5, we define the shadow state ip s N with iV G N as 



Cf 3 + C 2 5 a (-M for !/ = 4nz/ + 5 with n G [-JV, iV] 

elsewhere . 

(51) 

Clearly the shadow state is an almost periodic function and therefore f/>^r G $. Notice 
also that as iV — > oo, ^ approaches ^ G $*. The statement that the generic solutions 
have zero norm corresponds to the statement that 

(P^n) I p {^ 5 N))phys = (52) 

for 5 7^ <5 C where again 5 C = 4z/ [1 — mod(^^, 1)] according to (|50|). 

From the form of equation (}4T?j) it is clear that P 2 = P and hence the kinematical inner 
product can be used as the physical inner product. We can see explicitly that the solution 
1)50)1 is normalizable by calculating it's norm 

(E | E ) = ^2 p^Pu 

V 

CO 

= E W Gi) =1 < 53 > 



using a known formula for the sum over Bessel functions of integer order. 

The physical solution (|5Uj) is plotted in figure (JTJ). We see that despite the discrete 
nature of the solution, it approximates a continuous one which for large values of v would 
approximate the Wheeler-DeWitt equation. For positive values of v the solutions decays 
rapidly which would be expected for the Riemannian model which has no classically allowed 
region for an effective 5 positive cosmological constant. For negative v the solution oscillates 
with behavior indicating an effective negative cosmological constant. In the discussion that 
follows we will indicate why this occurs and the interpretation for the Lorentzian regime. 



5 Notice that the effect of the cosmological constant term depends on the orientation of the triad which is 
dynamical for our action. 
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FIG. 1: Physical solution (|5U|) . The dots indicate the values of the solution for where 4^ + 2 a^ 3 ^ ^' 
The actual solution would be zero elsewhere, however the solid lines are interpolated values designed 
to bring out the behavior of the solution more clearly. 

Furthermore, we can see that the classical singularity {y — 0) does not present a barrier to 
the evolution as the wave function can be evolved through it. 

V. DISCUSSION 

The simplicity of the Hamiltonian constraint used here has allowed for explicit calcula- 
tions to be performed. A rigorous definition of the physical scalar product was given and 
used to factor out spurious solutions to the constraint equations with zero physical norm. 
We wish to compare the results of the model with those of standard loop quantum cosmol- 
ogy. To do so, we observe that there is a simple relationship between the Riemannian and 
Lorentzian theories for homogeneous and isotropic models. Using this simple relation we can 
compare the two Hamiltonian constraints and the physical implications of the differences. 
In addition, we will discuss the cosmological predictions of such a model. 

The Lorentzian model and relation to standard LQC 

A Lorentzian model can be constructed if one uses general properties of the isotropic 
reduction of general relativity. With homogeneity and isotropy the curvature part of the 
Hamiltonian constraint changes sign while the cosmological constant remains the same when 
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going from the Riemannian sector to the Lorentzian one. Using this property the Hamil- 
tonian constraint for the Lorentzian constraint becomes Hl = —6A 2 + 3AE as opposed 
to the Riemannian constraint He = QA 2 + 3AE. Equivalently the same effect is obtained 
by changing the sign of the triad E (which amounts to sending v to — v in the solutions 
presented here). 

If we extend our model in such a fashion the difference equation becomes 

ijju+ivo - 2ip u + ip u - iuo = —^fvlll Aisipv (54) 

where we have included the various constants needed to compare to standard loop quantum 
cosmology. Here 7 is the Barbero-Immirzi parameter, l p = \fnh is the Planck length, and 
k = 8ttG. The difference equation for standard loop quantum cosmology is given by 

= A s 9 n (") V v (55) 

where the volume eigenvalues are V v = (|z/|7/ p /6) 3//2 , and sgn stands for the sign function. 
That the two constraints are asymptotically the same can be shown by using the large volume 
expansion for the volume difference coefficients V u+Uo — V v - Vo ~ (^ylp/Q) 3 ^ 2 3y/ui , o. Plugging 
this approximation into (|53|) then gives the difference equation ()54j) . The two models would 
thus share the same large volume semi-classical limit. 

The differences between the two difference equations can be traced back to the classical 
Hamiltonian constraints used. The constraint of standard loop quantum cosmology is 

H LQC = -— 2 A 2 sgn(E)^\E\ + -sgn(E)(\E\f 2 (56) 

whereas the constraint used here would be 

H = - — A 2 + -E. (57) 

It is readily seen that for positive values of the triad the two constraints are proportional 
to each other and thus are classically equivalent on the constraint surface where H = 0. 
The main qualitative difference between the two constraints is that they are not equivalent 
for negative orientations of the triad. This is indicated in the difference equation by the 
presence of the sgn{v) term on the right hand side of equation ()55|). The factor sgn(u) is 
put in by hand in standard LQC such that the constraint is classically symmetric for both 
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orientations of the triad. However, there is a freedom in extending the classical phase space 
to include negative orientations of the triad. We have shown that the action does not 
lead to the sgn(v) term. In the next subsection we discuss in more detail the implications 
of such a change. 

The other crucial difference between the constraints is that our constraint is self-adjoint. 
This allows us to carry out the group averaging technique to construct the physical inner 
product. The constraint of standard loop quantum cosmology is not, although self-adjoint 
constraints have been proposed j3, 3|- 



The triad orientation ambiguity 

We have stated the main qualitative difference between the simpler model lies in the clas- 
sical extension to negative orientations of the triad. This ambiguity, which is manifested by 
the presence of the sgn(u) term in the difference equation (jHHj) . has important consequences 
in the quantum theory. We now show that absence of the sgn(u) term — forced on us by 
the action (£[} — leads to many attractive features in the quantum theory. 

The physical wave function solutions of the two difference equations are shown in figure 
121 For large positive values of the triad (positive v) the solutions behave the same. However, 
for negative v the solutions are entirely different owing to the sgn{v) term. Because of that 
term, the solution in figure |2(b)| is symmetric and thus both orientations of the triad have 
the same semi-classical limit. On the other hand, the solution of the simplified constraint 
clearly does not have the right semi-classical limit for negative orientations of the triad. At 
first one might think that this is undesirable in the quantum theory since a large portion of 
the quantum configuration space does not have the right semi-classical limit. However, we 
see that the solution for negative triad is rapidly suppressed indicating a classically forbidden 
region. Thus, the universe would have very low probability for being in the region and would 
instead more likely be found in the semi-classical regime of positive v. Notice that because 
P 2 = P, the standard notion of probability amplitude can be associated to the physical 
wave function for the single partial observable E. Physically the cosmological constant term 
in (j57jl changes its sign when E goes to —E. The regions v > and v < correspond to an 
effective positive and negative cosmological constant respectively. 

It is precisely because of the absence of the sgn(v) term that the model can be solved 
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FIG. 2: Difference equation solutions for both the simplified constraint (a) and standard loop 
quantum cosmology (b). The solution of (a) is the physical solution. The solution in (b) is the 
one picked out by the dynamical initial condition. The solid line are the interpolated values. The 
actual physical solution of (a) would be zero except for the one sector indicated by the dots. For 
(b) the dynamical initial condition specifies that the wave function in the other sectors is the one 
that smoothly interpolates between the dotted solution. 

exactly. If we were to use a constraint with the sgn(v) term the path integral would not be 
solvable explicitly and it can be shown that there do not exist any kinematically normalizable 
solutions. The fact that the physical solution is kinematically normalizable simplified the 
construction of Ti. p h ys - Without kinematically normalizable solutions the path integral needs 
to be solved in order to construct the physical inner product. The normalizability of the 
solutions can be better understand by looking at the large volume behavior of the solutions. 
As in the Riemannian model the general solution of the simplified model is 



The Bessel Y solutions all diverge for negative v and the Bessel J solutions all diverge for 
positive v except for the one sector where the Bessel function order is an integer as it was 
in the Riemannian solution (j5()j) . The Bessel J solutions are plotted for large v in figure 
It is clear that only in the special sector 5 = 5 C does the solution not diverge. If we 
were to include the sgn(v) term then the oscillatory solutions would exist for positive and 




(58) 
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negative u, however we could not match a convergent solution for positive v with one for 
negative v. We thus see that without the sgniv) term, the model has very attractive features 
and the construction of the physical Hilbert space is greatly simplified since there exists a 
normalizable solution. 




3000 3200 3400 3600 3800 4000 3000 3200 3400 3600 3800 4000 
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(a) -i//(4i*) + 9/(2 7 3 /2A^ 3 ) S Z (b) -v/(lv Q ) + 9/(2^4) £ Z 

FIG. 3: Bessel J solution of ()58|) for two different sectors. The solution of (a) is for the sector 
where the Bessel function order is an integer. This solution converges both for large negative and 
positive v. The solution of (b) is for a different sector and diverges for large v. The solid lines are 
interpolated values. 

A question remains as to the relevance of the simplified constraint to that of standard 
LQC: are these special results specific to the simplified model? If we remove the sgn{y) 
term from the difference equation of standard LQC, it can be shown that the model with 
a cosmological constant does indeed exhibit similar behavior. To show this, the difference 
equation must be solved numerically since it is more complicated. Consider that we start 
at large negative v with some initial conditions of the wave function and we evolve forward 
toward the classical singularity. For large negative v the two independent solutions of the 
difference equation behave in a characteristic manner. The difference equation for large \v\ 
is approximated by the difference equation of our simplified model. Therefore, for large and 
negative v the two independent solutions behave as f+(v) ~ i/-( u+1 / 2 ) and f-(v) ~ v v ~ x l 2 
which can be shown from the asymptotic expressions of the solutions in equation (|58|). For 
generic initial conditions the solution will be a linear combination of /+ and /_. As we 
evolve toward the singularity /+ rapidly dominates over /_ . In this manner numerically the 
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contribution of the decaying solution /_ becomes negligible in comparison with the contri- 
bution of /+ thus effectively selecting the latter. We can perform the same trick starting 
at large positive v — where the independent solutions behave as g+(v) ~ (— l) 3 ^ iA -1 / 2 ) 
and g~(jy) ~ (— l) 3 ^ z/~(" +1 / 2 ) — and evolving backward to select the g + component. The 
two solutions can then be tested to see if they match somewhere in the region where the 
behavior is oscillatory. If the two can be matched then we have found a solution that decays 
both for positive and negative v. 

When this analysis is performed for the difference equation of LQC with a cosmological 
constant (without the sgn(v) term) the result is a match only in one sector which corresponds 
approximately to the one picked out by our simplified model and given by — u/(4vq) + 
9/(27 3 Z 2 Az/q) G Z. In the other sectors the matching cannot be achieved which implies for 
instance that a solution that decays toward — oo will evolve into one that diverges toward 
+00. Since the falloff behavior at large positive and negative v is sufficiently fast we can 
say that the solution that does not diverge is kinematically normalizable as in our simplified 
model. The solution is plotted in figure (0J). It has the same qualitative behavior as the 
normalizable solution of the simplified constraint (compare to figures p ( a) | and |3 ( a) \ . 

We are currently investigating whether a normalizable solution exists for a massless scalar 
field in the closed model. A recent paper has pointed out that in this model all the solutions 

n 

diverge either for large positive or negative v [3CJ. The question is raised that this implies 
that the wave function then predicts a large 'probability' for the universe being in those 
regions which are classically forbidden for this model. Notice however that any question 
about probabilities can only be addressed if a notion of physical inner product is provided. 
The physical inner product can solve this apparent problem by providing a non trivial 
probability measure that makes the physical states normalizable. We have shown for our 
model that most solutions have divergent behavior yet the physical inner product singles out 
a solution — all the other ill-behaved solutions have zero norm — which is finite everywhere 
and even kinematically normalizable. Therefore it is interesting to explore whether this 
might occur with the model of J^. If a normalizable solution can be found which decays 
for large volume then the problems raised would be shown to be nonexistent. Preliminary 
results indicate that one might be able to find normalizable solutions for this model again 
provided the sgn{v) term is removed. 
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FIG. 4: Normalizable solution to the difference equation of LQC (|55|) without the sgn(u) on the 
r.h.s. Here matter is only in the form of a cosmological constant. The solution is valued only in 
a single sector which approximately is given by the same sector of the simplified constraint where 
-*VOo) + 9/(2 7 3 /2A^) G Z. 

Cosmological implications 

We now discuss the cosmological implications of the Lorentzian model. A main question 
is whether or not the quantum theory cures the classical singularity. A classical singularity 
would exhibit itself in the quantum theory in two ways: a place where the difference equation 
breaks down or a place where the operator corresponding to the inverse scale factor is 
singular. For both these criteria the model presented here does not exhibit singular behavior. 
We see from the figure |2(a)| that the quantum evolution proceeds smoothly through the 
classical singularity iy = 0). From the standpoint of the difference equation the singularity 
point is not special. In addition the operator corresponding to the inverse scale factor is 
a kinematical operator which has been shown to be bounded in standard loop quantum 
cosmology jjilj]. Since the model here shares the same kinematical Hilbert space, the results 
will be the same. We furthermore see that the physical solution is valued only on discrete 
values of v corresponding to the special sector 5 C . The special sector depends on the value of 
the cosmological constant and only specific values will select the sector 5 = which passes 
through v = 0. Thus, it is most likely that the physical solution completely avoids the 
singularity. The singularity is thus cured in the model. 

If the model does not exhibit singular behavior then another question is what happens 
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when the universe approaches the singularity. The plot in figure 2(a) shows that the solution 
decays rapidly for negative values of v while it oscillates for positive values. Thus it is natural 
to interpret the region of negative v as a classically forbidden region and the singularity as 
a barrier off which the universe bounces. This is in contrast to standard loop quantum 
cosmology where semi-classical universes exist for both negative and positive v thus opening 
the possibility for the universe to tunnel through the singularity from one region to the 
other. Again, the main difference arises from the sgniv) term in the difference equation. 
From this we can make the interpretation that a classical collapsing universe approaching 
the singularity completely avoids it and bounces leading to an expanding universe. 

Finally, there is the issue of boundary proposals. In Wheeler-DeWitt quantization there 
exist two independent solutions for de Sitter space and ad-hoc boundary proposals are added 
to pick out a solution. In standard loop quantum cosmology the dynamical initial condition 
is used to pick out a solution. However, the dynamical initial condition only singles out 
a solution for the sector that passes through the singularity. Semi-classical arguments are 
then used to fix the wave function in the other sectors (see figure |2(bj| . In contrast we see 
that for the simplified model here, it is the physical inner product that selects a unique wave 
function. There is no need to supplement the quantum theory with an ad-hoc boundary 
proposal. This is not surprising for this model as the constraint equation is first order in 
the connection representation. The question remains as to whether this is relevant in more 
general situations. 



Further implications for LQC 

At this stage and in the context of the model presented here one question seems natural: 
what have we gained by setting up the theory on the non separable kinematical Hilbert 
space given by the Bohr compactification of the real line? The answer is clear in that 
considering quasi-periodic functions of A is the analog of considering cylindrical functions of 
the connection defined on arbitrary graphs in the full theory. Had we defined the quantum 
theory by using functions of a fixed periodicity (allowing v = uqu for n £ Z) we would 
have missed the physical solution Eq(A) and the physical Hilbert space would have been 
O-dimensional. Only in the special case when (3/2)A _1 z/^ 3 is an integer could we have 
found the physical state by starting with a formulation on a fixed 'lattice', i.e., where 
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v = uqti (see Equation (j5Uj)). In any other case using a fixed lattice would have resulted in 
a zero dimensional physical Hilbert space. In our model, the physical inner product selects 
a given set of 'graphs' by selecting a periodicity of the relevant modes. Thus the Bohr 
compactification is necessary to capture the correct physics. 

Another related issue is that the physical Hilbert space separates into orthogonal sub- 
spaces. In the standard quantizations of the Hamiltonian constraint in the literature the 
curvature term is written in terms of holonomies along paths whose length is determined by 
the parameter v . Due to this fact one introduces an intrinsic periodicity to the constraint. 
It is easy to see that the generalized projection operator associated to such a quantization 
satisfies the following property 

(v" \P\u')= JdT (u" | e- iAT | v') = (59) 

for v" — p 1 (fi Ai/qu for n G Z. In this manner the physical Hilbert space is separated into 
isolated sectors with no quantum interference between them. 
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FIG. 5: Wave function of standard LQC (dots) compared with the corresponding Wheeler-DeWitt 
solution (solid line) for positive cosmological constant. The deviation occurs when the extrinsic 
curvature is of the order of 7r/Vo which in this case occurs for large volumes. 

A direct consequence of this is the fact that in a given sector the physical solutions have 
a built in periodicity in that Au = 4z/ which implies that physically A < A max = 7r/z/ . 
This is a puzzling feature as A is directly related to the extrinsic curvature which classically 
does not have such a bound. In the model with a cosmological constant, classically the 
extrinsic curvature grows as the volume of the universe increases until at some volume it 
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reaches a critical value corresponding to A = A max = n/u . At values of v on the order 
of this critical value the behavior of the wave function changes dramatically and deviates 
from that expected from standard Wheeler-DeWitt quantization. This is a generic property 
that applies to all models of LQC (see figure EJ). However, it is in this region that both 
prescriptions should coincide. By decreasing the value of vq one can extend the region where 
LQC and Wheeler-DeWitt quantization agree, yet this par ameter is argued to be fixed by the 
full theory by considering the smallest area eigenvalue 2(| . It is surprising that a parameter 
arising from the fundamental discreetness at Planck scale should have such a important 
effect at large scales where the universe is expected to behave classically. 
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